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ABSTRACT 

When integrals in the standard Tremaine- Weinberg method are evaluated for the 
case of a realistic model of a doubly barred galaxy, their modifications introduced by 
the second rotating pattern are in accord with what can be derived from a simple 
extension of that method, based on separation of tracer's density. This extension 
yields a qualitative argument that discriminates between prograde and retrograde 
inner bars. However, the estimate of the value of inner bar's pattern speed requires 
further assumptions. When this extension of the Tremaine- Weinberg method is applied 
to the recent observation of the doubly barred galaxy NGC 2950, it indicates that the 
inner bar there is counter-rotating, possibly with the pattern speed of —140 ± 50 km 
s _1 arcsec -1 . The occurrence of counter-rotating inner bars can constrain theories of 
galaxy formation. 
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1 INTRODUCTION 

Bars within bars appear to be a common phenomenon in 
galaxies. Recent surveys indicate that up to 30% of early- 
type barred galaxies contain such double bars (Erwin & 
Sparke 2002; Laine et al. 2002). Inner bars remain distinct 
in near infrared (Wozniak et al. 1995), therefore fairly old 
stars must contribute to their light. The relative orientation 
of the two bars in doubly barred galaxies is random, there- 
fore it is likely that the bars rotate with different pattern 
speeds. 

The origin of multiply barred systems remains unclear. 
A bar takes away angular momentum from gas very effi- 
ciently, and in the young Universe bars might have been 
responsible for the early rapid growth of the massive black 
holes (Begelman, Volenteri & Rees 2006). If the innermost 
parts of galactic discs formed first, then early instabilities 
there might have lead to the formation of small-scale bars, 
which may be surviving in the present-day Universe as nu- 
clear bars, often nested inside larger, outer bars, that might 
have formed later. Studying the dynamics of nested bars can 
therefore help us to understand the formation of galaxies and 
of their central massive black holes. 

Our understanding of how double bars are sustained 
has significantly improved in the recent years. Maciejewski 
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& Sparke (1997, 2000) have developed a formalism that en- 
ables finding families of stable regular orbits in such systems. 
Stable regular orbits are robust structures, which define the 
shape of the galaxy, and therefore they can serve as a back- 
bone for double bars. So far they have been analyzed for the 
case when the pattern speed of the inner bar is higher than 
that of the outer bar (Maciejewski & Sparke 2000). Using 
N-body simulations, Rautiainen et al. (2002) confirmed that 
stars get trapped around these orbits, and form long-lasting 
doubly barred systems. However, in this scenario the kine- 
matics of the inner bar is not a scaled-down copy of that of 
the main bar, since the inner bar cannot extend to its corota- 
tion. Other numerical simulations have shown that systems 
of two counter-rotating bars are also possible (Sellwood & 
Merritt 1994; Friedli 1996) , and systems with secondary bars 
rotating slower than the outer, main bars, have never been 
excluded on theoretical grounds. 

Various dynamical scenarios for doubly barred galax- 
ies may have different implications for the evolution of the 
galactic centres. In order to discriminate between them, one 
should measure the pattern speed of the inner bar. If only 
one pattern speed is present in the system, the Tremaine- 
Weinberg (1984) method allows to derive it using a set of 
simple kinematical measurements. Recent kinematical ob- 
servations of the doubly barred galaxy NGC 2950 (Corsini, 
Debattista & Aguerri 2003; hereafter CDA03) are inconsis- 
tent with one pattern speed there. The observations are sug- 
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gestive of another pattern speed in the area of the inner bar. 
CDA03 attempted to estimate it, but they concluded a wide 
range of values, consistent with a fast-rotating prograde sec- 
ondary bar as well as with a retrograde one. 

In this paper, I show that a simple extension of the 
Tremaine- Weinberg method to multiple pattern speeds, 
based on the separation of tracer's density, is sufficient to 
discriminate between prograde and retrograde inner bars. 
Similar extension has been already considered by CDA03, 
but they did not explore its implications. In Section 2, I out- 
line the extended method, and I show that it immediately 
gives a qualitative information about the sign of rotation of 
the inner bar. In Section 3, I calculate the integrals in the 
original Tremaine- Weinberg method for the case of a real- 
istic doubly barred galaxy, and I show that their deviation 
from values for a single rotating pattern is the same as pre- 
dicted by the extension to the Tremaine- Weinberg method 
proposed in Section 2. In Section 4, I examine one method 
to measure the pattern speed of the inner bar, and I show 
that it recovers the pattern speed in the model with an 
acceptable accuracy. Limitations of the extended method, 
other attempts to estimate multiple pattern speeds using the 
Tremaine- Weinberg formalism, and consequences of counter- 
rotating inner bars for the evolution of galactic centres are 
discussed in Section 5. 

When the extended Tremaine- Weinberg method pro- 
posed here is applied to NGC 2950, it indicates that the 
pattern speed of the inner bar there is smaller than the pos- 
itive pattern speed of the outer bar, and that it is most likely 
negative, i.e. the inner bar is counter-rotating with respect 
to the outer bar and to the disc (see also Maciejewski 2004) . 



2 A SIMPLE EXTENSION OF THE 

TREMAINE- WEINBERG METHOD TO 
MULTIPLE PATTERN SPEEDS 

The Tremaine- Weinberg method is designed for one pattern 
speed, which it derives from the luminosity centroid and the 
luminosity-weighted line-of-sight velocity of a chosen tracer 
moving in the galaxy's potential. The method rests on three 
assumptions: the disc of the galaxy is flat, it has a well- 
defined pattern speed, and the tracer obeys the continuity 
equation. Under these assumptions, the surface density of 
the tracer, T,(x,y,t), can be written as 

V(x,y,t) = £(R,<p-n P t) (1) 

where flp is the pattern speed, (x,y) are Cartesian coordi- 
nates in the disc plane, (R, ip) are polar coordinates there, 
centred on the galactic centre, and t is time. 

When another pattern speed is introduced, the two pat- 
terns cannot rotate rigidly one through another (Louis & 
Gerhardt 1988; Sridhar 1989; Maciejewski & Sparke 2000), 
and in principle one cannot split the right-hand side of (1) 
into components with different pattern speeds. However, if 
a rough separation of patterns is possible, and if one ne- 
glects secular evolution, this system is periodic with period 
P = tt/(Qs — Hs), where fis and Q, B are pattern speeds of 
the two bars. Consequently, the surface density of the tracer 
can generally be written as 

T.{x,y,t) = £ B (R,<p - Sl B t,t\P) + t s (R,f -Q s t,t\P), (2) 



where t\P denotes dependence on time with periodicity P. In 
Appendix A, I evaluate correction to the Tremaine- Weinberg 
integrals that arises because of this periodic oscillation of 
realistic double bars. This correction turns out to be small, 
and I will neglect it in the following argument. 

Once periodic oscillations of nested bars are neglected, 
(2) gets reduced to 

X(x,y,t) = £ B (R,<p - n B t) + £ s (R,<fi - n s t). (3) 

Twofold interpretation of such a separation of tracer's den- 
sity is possible. Either the tracer in the disc (e.g. stars) can 
be divided into subgroups belonging to two unchanging pat- 
terns, each rotating with a constant pattern speed, or one 
can separate radial zones in the galactic disc that rotate with 
constant pattern speeds, with no net exchange of the tracer 
between the zones. Neither of these interpretations is con- 
sistent with the dynamics of galaxies, because patterns ro- 
tating one through another change with time, and tracers of 
each pattern are often present at the same radius. However, 
as I show in Appendix A, these inconsistencies are likely to 
be small, and the separation (3) can still be approximately 
valid. 

Below I will show that although the extension of the 
Tremaine- Weinberg method based on the separation of 
tracer's density (3) cannot recover the two pattern speeds, it 
can yield a qualitative prediction of how the integrals in the 
standard Tremaine- Weinberg method change when a second 
rotating pattern is introduced. In Section 3, I will show that 
this change, calculated properly for a realistic model of dou- 
bly barred galaxy, with bars oscillating in time, and with 
tracers of each bar overlapping, is the same as predicted by 
the extended method proposed here. Thus the separation 
of tracer's density may help in interpreting the Tremaine- 
Weinberg integrals in the models and in the observed galax- 
ies. 

Once tracer's density is separated according to (3) , one 
can write two separate continuity equations. Further deriva- 
tion is identical to the one performed by Tremaine & Wein- 
berg (1984), and it can be conducted for each tracer's density 
separately, leading to two equations: 

/ + oo 
E B (X,Y 3llt )XdX = 
-OO 

/ + OG 
X B (X,Y sHt )V B os (X,Y 3Ut )dX, (4) 
-OO 

/ + oo 
T, s (X,Y slit )XdX = 
-oo 

/ + CO 
Xs(X,Y sUt )Vt s (X,Y sUt )dX, (5) 
-OO 

where (X,Y a m) are coordinates on the sky (X running par- 
allel to the line of nodes and Y s i it being the offset, per- 
pendicular to the line of nodes, of the slit along which the 
integration is performed), V los is the observed line-of-sight 
velocity, and i is the inclination of the disc. Let the tracer 
be made of stellar-light emission; then one can define for 
each bar component, as well as for the whole galaxy, the 
luminosity-density 

/ + oo 
£o(X, Y s ut)dX, 
-OO 
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Figure 1. Schematic drawing of the positions of the bars ob- 
served in NGC 2950 and of an exemplary location of the slit. 
Two solid ellipses outline the bars, and the dotted ellipse marks 
the boundary used in the derivation in Section 4. 

the luminosity centroid 
1 [ +x 

Xo = — / Eo(X,Y M t)XdX, 

•/ — oo 

and the luminosity-weighted line-of-sight velocity 
1 f + °° 

Vo = — Eo(X, Y Mt ) Vi os (X, Y sllt )dX, 

where the index o is B for the outer bar, S for the inner 
bar, and tot for the whole galaxy. Obviously, only values 
with index tot can be observed, for example by placing a slit 
parallel to the line of nodes (Fig.l), and deriving along it 
the values of luminosity, centroid and line-of-sight velocity, 
where all integrals include integration over the width of the 
slit. With the definitions above, the sum of (4) and (5) takes 
a form 

Q B sin i F B X B + Qs sin i F S X S = F B V B + F s Vs, (6) 

where F<> = L /Lt o t- CDA03 gave a similar equation, but in 
their convention there are no fractions F . Note that equa- 
tion (6) is not sufficient to derive two pattern speeds Qb,s by 
measuring the luminosity centroid Xtot and the luminosity- 
weighted line-of-sight velocity Vtot along the slits. Although 
FbVb + FsVs = Vtot, and it can be measured directly, the 
pattern speeds are combined with the unknown Fb,s and 
Xb,s- 

However, the form of equation JJj), together with the 
morphology of the doubly barred galaxy, whose pattern 
speeds we want to measure, already reveals some informa- 
tion about the relation between the two pattern speeds. The 
line of nodes passing through the centre of the galaxy sets a 
division of this galaxy into four quadrants. Consider for ex- 
ample the case when the two bars lie in opposite quadrants 
of the galaxy (Fig.l). This is the case of NGC 2950. Even if 
we cannot measure Xb or Xs, we know that they are always 
of opposite signs, no matter how the slit is placed (parallel 
to the line of nodes). Then if Qb and Qs are of the same 
sign, adding the inner bar should bring the sum on the left of 
(jSJ closer to zero (or even through zero for strong fast inner 
bars) when compared to the contribution of the outer bar 



alone. This implies that the observed \Vtot \ should be smaller 
when entering the region of the inner bar than that interpo- 
lated from the outer bar. To the contrary, CDA03 observe 
\Vtot\ increasing around this transition region in NGC 2950, 
consistent with Qs having the sign opposite to Qb, i.e. with 
the counter-rotating secondary bar. 

2.1 Immediate condition for the sign of Qs Qb 

In the argument above, no reference point has been fixed, 
since the sums on each side of l|SJ get modified by the in- 
troduction of the secondary bar. We do not know the con- 
tribution from the outer bar in the region where the inner 
bar is present, and interpolation can be misleading. Here I 
show that © allows to tell whether Qs > Qb or Qs < Qb- 
Equation JSJ can be rewritten as 

Q B Xt t+X s Fs(Qs - Qb) = V to t/smi. (7) 

This form shows how the relation between the integrals Vtot 
and Xtot gets modified by the presence of the secondary bar: 
the second term in the sum is a correction term. For slits 
outside the inner bar this correction is zero (because Es is 
zero there), and for those slits the observed values of Xtot 
and Vtot should lie on a straight line of inclination Qb sini. 
When the slit passes through the inner bar, Xs 7^ 0. If 
for such slits \Vtot\ is larger than the values given by the 
linear relation (7) without the correction factor, then the 
correction has to have the same sign as QsXtot- For the bars 
that lie in opposite quadrants, like in NGC 2950 (Fig.l), Xb 
and Xs are of opposite signs, and Xtot has the same sign as 
Xb throughout the galaxy (fig. 3 in CDA03). Therefore the 
two components of the sum in (7) can only be of the same 
sign when Qb has opposite sign to Qs — Qb- If we take a 
convention that Qb > then Qs < Qb- 

This argument can be extended to the case when Xtot 
changes sign for slits passing through the secondary bar. In 
NGC 2950, the values of Xtot for these slits are very small, 
which is consistent with the observed geometry of bars. A 
similar argument can be applied to galaxies with bars in the 
same quadrants, where Xtot for innermost slits is not that 
small. In this case, Xb, Xs and Xtot are all of the same sign. 
If an increase of \Vtot\ is observed in the region of the inner 
bar, Qb must now have the same sign as Qs — Qb, which 
means Qs > Qb for the assumed Qb > 0. Thus we see that 
the same increase of \Vtot\ in the region of the inner bar can 
either indicate the inner bar rotating slower or faster than 
the outer bar, depending on the relative orientation of the 
bars. 

2.2 Condition for counter-rotation 

In some cases, an argument can be made about corotation 
or counter-rotation of the inner bar with respect to the in- 
ertial frame. It relies on the variation of the integrals in 
the Tremaine- Weinberg method with the distance of the slit 
from the centre of the galaxy, Y s i it . This argument can be 
made under an assumption that \XbFb\ does not increase 
when we march with the slit through the galaxy toward the 
line of nodes (in the direction marked by an arrow on the 
left of Fig.l). This is a justified assumption, since Fb de- 
creases inward for the very reason of the introduction of 
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Figure 2. Loops (solid lines) in a realistic model of doubly barred galaxy (Model 2 of Maciejewski &; Sparke 2000) in the linear 
approximation. Dashed ellipses outline the bars: the density of each bar drops there to 0.1 of its central value. Dotted lines mark the 
direction along which 'slits' are placed, for which the integrated values of Xtot and Vtot are derived. The bars are aligned in the left 
panel, and at 57r/8 rad angle in the right panel. Units on axes are in kpc. 



the secondary bar, and \Xb \ should not go up inward, since 
early-type galaxies like NGC 2950 have 'flat' bars with a 
nearly constant surface brightness as a function of radius 
(Elmegreen et al. 1996). 

Consider again equation ©. Normally, for slits that 
avoid the secondary bar, \Vtot\ decreases inward (i.e. when 
shifting the slit in the direction given by the arrow in Fig.l). 
This is also the case in NGC 2950 (fig.3 in CDA03). If then 
\Vtot\ increases inward when the slit reaches the secondary 
bar (like in NGC 2950), this can be caused by either of 
the two components of the sum on the left of I already 
argued that the first component cannot be the cause, since 
Fb decreases, \Xb \ is unlikely to increase, and Qb = const. 
The second component can only cause the increase of the 
sum when it is of the same sign as the first one. But again, 
Xb and Xs are of opposite sign. Therefore Os and Qb have 
to be of opposite sign, too. 

One can repeat this reasoning for a galaxy with bars in 
the same quadrants, to show that in that case | Vtot | increas- 
ing in the region of the inner bar indicates that the inner 
bar is corotating in the inertial frame. Similar to Section 2.1, 
the conclusion about the sense of rotation of the inner bar 
depends on the relative orientation of the bars. 



3 INTEGRALS IN THE 

TREMAINE- WEINBERG METHOD 
CALCULATED FOR A REALISTIC MODEL 
OF A DOUBLY BARRED GALAXY 

Regular motion of a particle in a potential of a doubly 
barred galaxy has two frequencies associated with it, each 
related to one of the bars, in addition to the frequency of 
its free oscillations (Maciejewski & Sparke 1997). In the lin- 
ear approximation this motion corresponds to epicyclic os- 



cillations with these frequencies around the guiding radius. 
Particles with the same guiding radii are bound to closed 
curves (loops) that oscillate in the pulsating potential of 
double bars. Loops are also observed in nonlinear analysis 
(Maciejewski & Sparke 2000; Maciejewski & Athanassoula 
2006). 

The orbital approach directly indicates that the motion 
of a particle associated with one bar has always a compo- 
nent coming from the other bar. Amplitudes of the oscil- 
lations can be easily evaluated in the linear approximation 
(Maciejewski 2003), giving particle's position and velocity at 
each relative position of the bars. In Fig. 2, I plot the loops 
populated by particles moving in the potential of a doubly 
barred galaxy constructed by Maciejewski & Sparke (2000; 
Model 2). This is a realistic potential, since it admits or- 
bits that support the outer bar, as well as orbits supporting 
the inner bar, throughout the extent of each bar. The two 
bars in this model rotate in the same direction with pattern 
speeds Qb = 35 km s _1 kpc -1 and Qs = 110 km s _1 kpc -1 . 
I plot the loops for two relative orientations of the bars: the 
bars parallel (Fig. 2, left), and the bars at the angle of 57r/8 
rad, the value similar to that observed in NGC 2950 after 
deprojection (Fig. 2, right). The loops change shapes as the 
bars rotate one with respect to another, and loops associated 
with one bar intersect the loops associated with the other 
bar. Thus formally one cannot perform here the separation 
of tracer's density postulated in (3). 

However, changes in the shapes of the bars, and zones 
where both tracers coexist, may be small, and then, what 
formally prohibits the separation, may turn into a higher- 
order correction to it (see also Appendix A). In order to 
check whether it is the case for the realistic Model 2 (Ma- 
ciejewski & Sparke 2000), I calculated the positions and ve- 
locities of some 10 6 points on 200 loops in that model, from 
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Figure 3. Values of the integrals in the Tremaine- Weinberg method, calculated from the epicyclic formulae for the positions and velocities 
of particles in the realistic model of doubly barred galaxy from Fig. 2. The dot-dashed line corresponds to the model without the inner 
bar, the dashed line is for the model with two bars aligned, as in the left panel of Fig. 2, while the solid line is for the model with two bars 
located as in the right panel of Fig. 2. This figure can be directly compared with fig.3 in CDA03. (a) The kinematic integrals Vtot, defined 
in Section 2, as a function of the slit offset Y s m with respect to the centre of the galaxy, (b) The photometric integrals Xtot, defined 
in Section 2, as a function of the slit offset Y sHt . (c) Vtot as a function of Xtot- (d) Pattern speed in the model, O pa ttern = Vt t/Xt t, 
derived separately for each position of the slit Y a ut- 



which I obtained the centroids X to t and the line-of-sight ve- 
locities Vtot along slits placed at the same relative angle to 
the outer bar as in the CDA03 observations of NGC 2950 
(dotted straight lines in Fig. 2). For clarity, only 28 loops are 
displayed in Fig. 2, out of 200 used in the calculations. The 
200 loops were populated with particles, so that a smooth 
density distribution in the bars is recovered. Several recipes 
have been adopted to populate loops with particles, in order 
to make sure that the results below do not depend on the 



way in which the model is constructed. In the model, I use 
deprojected velocities, which replaces Vtot/ sin i by Vtot, and 
I measure Y a m in the plane of the galaxy. 

In Fig. 3b, I plot the value of the centroid, X to t, as a 
function of the offset of the slit from the galaxy centre, Y 3 u t . 
When the inner bar is placed at the position similar to the 
one observed in NGC 2950, the value of X tot (solid line) is 
closer to zero in the region of the inner bar (|Y],iit| < 1), than 
when the inner bar is absent (dot-dashed line). However, 
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Xtot does not change the sign in the region of the inner 
bar, which is also observed in NGC 2950 (fig.3 in CDA03). 
This indicates that the light integrated along a slit passing 
through the inner bar is still dominated by the outer bar. 
In the case of the two bars parallel (Fig.2, left panel) there 
is no such decrease of X to t in the region of the inner bar 
(dashed line in Fig. 3b) . 

In Fig. 3a, I plot the luminosity-weighted line-of-sight 
velocity, Vtot, as a function of the same offset of the slit 
from the galaxy centre, Y s nt, as in Fig. 3b. When the in- 
ner bar, rotating in the same direction as the outer bar, but 
faster, is placed at the position similar to the one observed in 
NGC 2950, the value of Vtot (solid line) rapidly approaches 
zero in the region of the inner bar (|Y s j»t| < 1), and in most 
of this region it has the sign opposite to the value of Vtot for 
slits not passing through the inner bar (|Y s (jt| > 1). This be- 
haviour of Vtot is opposite to that observed in NGC 2950 by 
CDA03. For the model of galaxy with Qs > Qb, changes of 
Vtot with Y s ut observed by CDA03 can only be reproduced 
when both bars lie in the same quadrants of the galaxy. For 
the particular case of bars parallel, presented here, |Vtot| in- 
creases in the region of the inner bar (dashed line in Fig. 3a), 
relative to the case when there is no inner bar (dot-dashed 
line). This dependence of deviations in Vtot on the relative 
position of the two bars is exactly as expected from the 
simple extension of the Tremaine- Weinberg method derived 
Section 2, based on the separation of tracer's density (3). 
Namely, if the inner bar rotates faster than the outer bar 
then, when both bars are in the same quadrants (defined 
by the line of nodes and the centre of the galaxy), \Vtot\ 
increases in the slits passing through the inner bar, but it 
decreases, when the bars are in the opposite quadrants. 

The plot of Vtot as a function of Xtot is presented in 
Fig. 3c. One may attempt to fit a second straight line in the 
region of the inner bar, but this fit will not yield the actual 
pattern speed of the inner bar. In Fig. 3d, I plot the ratios of 
Vtot I Xtot as estimators of the derived Qpattem- For a single 
bar (dot-dashed line), there is one pattern speed indepen- 
dent of the offset of the slit, and equal to the one assumed 
in the model. If there are two pattern speeds in the system, 
the inner bar does not alter significantly the pattern speed of 
the outer bar derived from the slits that do not pass through 
this inner bar. However, the induced deviations are system- 
atic: pattern speed of the outer bar is spuriously increased 
when the bars lie in the same quadrants, and decreased when 
they lie in quadrants opposite. The effect is small though, 
below 5%. For the slits passing through the inner bar, the 
derived Q pa ttem rapidly decreases if the bars lie in opposite 
quadrants (solid line), and it becomes negative, not giving 
any information about the pattern speed of the inner bar. 
If the bars lie in the same quadrants (dashed line), the de- 
rived Qpattem is larger than Qb for the slits passing through 
the inner bar, but it varies with Y s m, and never reaches the 
value of fis = HO km s _1 kpc -1 assumed in the model, 
remaining at slightly over half of that value. 



4 A METHOD TO DERIVE THE VALUE OF Q s 

The main goal of this paper is to show that with a sim- 
ple extension of the Tremaine- Weinberg method to multi- 
ple pattern speeds one can derive at least rough qualitative 
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Figure 4. Model values of the integrals that determine the linear 
regression in (10) with the expected slope of Qg — Values of 
tot I-x ^tot(X, Y s ut)XdX are plotted against Vtot — VlsX-to 
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information about the secondary pattern rotation. For the 
data on NGC 2950 presented by CDA03 this information is 
that the inner bar rotates in the opposite direction than the 
outer bar. On the other hand, getting the numerical value of 
the pattern speed of the secondary bar, Qs, may not be pos- 
sible without making additional assumptions. Here I analyse 
one possible method to calculate Qs- 

If the tracer of the secondary bar does not extend be- 
yond some radius Ro in the galaxy plane (outlined in projec- 
tion by the dotted ellipse in Fig.l), then (5) can be rewritten 
as 

Q s smi ( ° T, s (X,Y Mt )XdX = 
J-x 

f + X 



Y,s(X,Y Mt )V s os (X,Y Mt )dX, 



(8) 



where R 
to 

Q B X to t + 



(Y sl 



t/cosif + Xq. 



Summing (4) and (8) leads 



Qs - Qe 



x 



T, s (X,Y Mt )XdX = ^SL. (9) 

Sill! 



However, aside for Qs, we still do not know the integral in 
(9). It can be approximated when one assumes that only 
tracers of the inner bar and of the axisymmetric component 
are present inside the dotted ellipse in Fig.l, i.e. that the 
outer bar is almost axisymmetric in this region. Note that 
this is a strong and poorly founded assumption. However, if 
we take it, and since the contribution of the axisymmetric 
component to the integral in (9) cancels out, one can sub- 
stitute there the observable E tot for the unknown Eg, and 
rewrite (9) as a linear regression of Qs — Qb- 

° T, to t{X, Y Mt )XdX = ^L- Q B X t0 f (10) 
sin i 

-x 



Qs - Qi 
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Figure 5. Values of f " Etot(X, Y s ut)XdX as a function 

l^tot J — Aq 

of R = y/(Y slit / cosi) 2 + X 2 for NGC 2950. The curves arc 
derived from the R-band image of NGC 2950 (Erwin & Sparkc 
2003), for slits placed parallel to the line of nodes, and offset from 
it by values given on the right of the plot in arcsec. Two dotted 
vertical lines indicate a plateau range, common for all the curves. 

Thus Qs — Qb can be obtained as a slope of a straight line 
fitted to the data from slits passing through the inner bar. 
A similar equation has been used by CDA03 to estimate Qs, 
but here the coefficient at Qs — Qb is defined differently. 

I tested equation (10) on the model examined in Section 
3, for the position angles of the bars like the ones observed in 
NGC 2950. As in Section 3, in the model I use deprojected 
velocities, which replaces Vtot/ sini by Vtot, and I measure 
Y s ut in the plane of the galaxy. On the basis of the shapes 
of the loops, I chose Ro = 1.1 kpc for this model, and I fixed 
Qb — 35 km s" 1 kpc -1 . In Fig. 4, I plot the data points 
for the regression (10). They follow a straight line well, ex- 
cept for the points around zero values. However, these points 
come from Y s u t ~ Ro, where the integration is over a small 
number of particles, hence random error is large there. If we 
exclude these points, the derived Qs — Qb oscillates between 
60 and 70 km s" 1 kpc -1 for |Y s j it | < 0.75 kpc. This gives 
Qs between 95 and 105 km s -1 kpc -1 , consistent with the 
input value of 110 km s" 1 kpc -1 , with an error of about 
10%. 

In order to apply the same method to the observed data 
on NGC 2950 from CDA03, the value of R can be esti- 
mated by checking how the integral in (10) changes with 
varying Ro — it should have a plateau in the transition 
area between the bars. I used the R-band image of NGC 
2950 (Erwin & Sparke 2003) to calculate this integral and 
Ltot for each of the slits passing through the secondary 
bar, placed at positions reported by CDA03. The values of 
ITT f-x ^tot(X, Y 3 ut)XdX as a function of Ro are plotted 
in Fig. 5. All curves indeed have a plateau around the same 
Ro of 5-6 arcsec. It is located just outside a nuclear stellar 
ring at ~ 4.2 arcsec, reported by Erwin & Sparkc (2003), 
which is likely circular in the plane of the galaxy. Thus the 
method proposed here may be particularly well suited for 
NGC 2950. 

Each curve in Fig.5 has a characteristic U-shape, whose 



depth depends on the relative contribution of the inner bar 
to the integral in (10), and therefore it vanishes outside the 
inner bar, and also very close to the line of nodes, where 
the integral tends to be zero. From the R-band image, the 
depth is largest for slits offset by ~ ±1 arcsec. The linear 
fit to (10), with the data for Vtot and X to t from each slit 
passing through the inner bar provided by CDA03, reveals 
a large negative value of Qs — Qb, about —150 ± 50 km 
s" 1 arcsec -1 . For Qb — 11 km s" 1 arcsec -1 , measured by 
CDA03, this gives Q s ~ -140 ± 50 km s" 1 arcsec" 1 . This 
result, although with much poorer grounds than the quali- 
tative arguments from Section 2, reaffirms counter-rotation 
of the inner bar. Note that the unrealistically large value of 
Qs derived with this method is likely an overestimate: in or- 
der to get it, the integral J X ^ o T,tot(X,Y s i it )XdX has been 

substituted for f*° T, s (X,Y Mt )XdX. The substituted in- 
tegral is most likely smaller than the original one, because 
the contributions of the two bars are of opposite signs. 



5 DISCUSSION 

The original Tremaine- Weinberg method, and its extension 
proposed here, are based on the continuity equation applied 
to a tracer moving in the gravitational field of the galaxy. 
If old stars are taken as the tracer, the continuity equation 
is well satisfied globally. If tracers associated with each bar 
are separated as proposed in (2), continuity of each tracer 
can only be violated when there is a net flux of stars from 
one bar to the other, consistent over many rotations of the 
bars. This corresponds to a secular strengthening of one bar 
at the cost of the other, hence such a system is no longer 
periodic. However, if the mass transfer is slow, the analy- 
sis presented in this paper is still applicable. This can be 
supported by an argument similar to that presented in Ap- 
pendix A, but for densities of the bars monotonically chang- 
ing. Moreover, studies of orbits in self-consistent models of 
double bars show that mass transfer between the bars is 
likely to be small, because orbits that are trapped and os- 
cillate around one or the other bar, populate large fraction 
of phase space (Maciejewski & Sparke 2000; Maciejewski & 
Athanassoula 2006). 

Note that the separation of tracers in (2) is only formal, 
and (2) applies to any system that is periodic with period 
P. However, in order to follow the Tremaine- Weinberg for- 
malism for each tracer, one has to approximate each tracer 
as a solid-body rotator, which leads to (3). The analysis of a 
realistic doubly barred galactic potential, presented in Sec- 
tion 3, indicates that although this approximation is incor- 
rect in a rigorous sense, deviations from solid-body rotation 
are small, and they only contribute to higher-order terms in 
the Tremaine- Weinberg integrals. Namely, even if each par- 
ticle in the system oscillates with frequencies related to both 
bars, for most of the particles one frequency dominates, and 
the separation can be performed. Although the presence of 
two forcing frequencies leads to noticeable pulsation of the 
inner bar (e.g. Rautiainen et al. 2002), the magnitude of ad- 
ditional terms in the Tremaine- Weinberg integrals, that this 
pulsation gives rise to, is much smaller than the magnitude 
of the leading terms, as argued in Appendix A. The radial 
separation of the tracers, necessary for the estimate of the 
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pattern speed of the inner bar in Section 4, can also be pos- 
sible, if the zones, where two tracers coexist, occupy small 
fraction of the galaxy. The larger this fraction, the poorer is 
the estimate of the inner bar's pattern speed. Further test- 
ing of this method will require a fully self-consistent model 
with an inner bar. 

This paper shows that the change in value of the inte- 
grals Xtot and Vtot in the original Tremaine- Weinberg for- 
malism, caused by the second pattern speed, depends on the 
relative position of the bars. In particular, I showed that fit- 
ting another straight line in the Xtot — Vtot diagram to the 
data from slits passing through the inner bar does not yield 
the pattern speed of this bar. Therefore recent derivations 
of multiple pattern speeds based on such analysis (Rand & 
Wallin 2004; Hernandez et al. 2005) have to be treated with 
caution. 

The retrograde rotation of the inner bar in NGC 2950 
proposed here is inconsistent with the first estimate of Q,s 
by CDA03. However, in that estimate it is assumed that the 
inner bar dominates the light in the slits passing through it. 
This cannot be correct, since X to t and Xs are of opposite 
sign there, which indicates that the light from the outer bar 
still dominates in the slits passing through the inner bar. 
The second estimate of CDA03, like the estimate given in 
this paper, indicates that the inner bar is counter-rotating. 

The obvious next step in verifying counter-rotation of 
the inner bar in NGC 2950 would be an examination of 
orbital structure of a galaxy with a retrograde inner bar, 
following the method presented in Section 3. However, a 
bar counter-rotating with respect to its stellar disc has to 
be built out of so called X4 orbits, which are normally very 
close to circular (see e.g. Sellwood & Wilkinson 1993). In 
terms of the linear approximation of Section 3, the ampli- 
tude of oscillation around the guiding radius generated by a 
retrograde bar is very small (see equations (12) and (13) in 
Maciejewski 2003). Thus a construction of a retrograde bar 
in a prograde stellar disc is unlikely. This is consistent with 
N-body simulations, in which multiple pattern speeds form 
naturally in stellar discs (Rautiainen et al. 2002) — these 
simulations do not recover retrograde pattern speeds. 

This leads to another possibility, namely that the inner 
bar is formed out of an inner disc that counter-rotates with 
respect to the outer disc. Numerical results indicate that 
two counter-rotating bars, formed in two counter-rotating 
stellar discs that overlap each other, can survive in galaxies 
for many rotation periods (Sellwood & Merritt 1994; Friedli 
1996). If this is the case, one should find out whether a 
retrograde inner disc is consistent with the rotation curve 
and velocity dispersion in NGC 2950. The data presented 
by CDA03 do not indicate counter-rotation in the innermost 
few arc-seconds of NGC 2950. 

If further kinematical studies of NGC 2950 confirm 
counter-rotation in its innermost parts, then extending the 
Tremaine- Weinberg method along the lines proposed here 
may yield an efficient way of detecting counter-rotation in 
galaxies with nested bars. Thus far all detections of counter- 
rotation in galaxies come directly from the observed veloc- 
ity fields obtained with the long-slit (e.g. Kuijken, Fisher & 
Merrifield 1996) or integral-field (e.g. Emsellem et al. 2004) 
spectroscopy. However, if the counter-rotating population is 
small, it may not be recognized with those methods. On the 
other hand, the Tremaine- Weinberg method goes beyond 



the raw observed velocity field by finding integrals that are 
useful in detecting a rotating pattern. In the presence of 
rotating patterns counter-rotation should be spotted more 
easily with this method. 

The use of this method, like of the original Tremaine- 
Weinberg method, is limited to galaxies with bars consid- 
erably inclined to both the major and minor axes of the 
disc. NGC 2950 fulfils this requirement particularly well, 
but there is a number of other galaxies to which this method 
can be readily applied. For example NGC 3368, NGC 3941, 
NGC 5365, NGC 5850, and NGC 6684 all have the po- 
sition angles of both bars separated by 20 to 70 degrees 
from the position angle of the line of nodes (see the cat- 
alogue by Erwin 2004 for the parameters). The bar inside 
the oval in NGC 3081 has also parameters favourable for 
this method. In addition, a modification of the Tremaine- 
Weinberg method along the lines proposed in Section 4 can 
be applied to the interiors of nuclear rings that often host 
nuclear bars (e.g. NGC 1097, NGC 6782). 

The fraction of counter-rotating inner bars can con- 
strain theories of galaxy formation and evolution. Currently 
the most accepted view on the origin of the inner bar is 
that it forms through instabilities in gas inflowing along the 
outer bar (Shlosman, Frank & Begelman 1989). However, if 
inner bars form on early stages of galaxy assembling, and 
outer bars form from material that settled on the galaxy 
later, the spins of the two bars may be unrelated, leading 
to a large fraction of counter-rotating inner bars. Determin- 
ing this fraction may help to distinguish between these two 
evolutionary scenarios. 



6 CONCLUSIONS 

In this paper I presented an argument that a simple exten- 
sion of the Tremaine- Weinberg method to multiple pattern 
speeds can provide an information about the sense of ro- 
tation of the inner bar in doubly barred galaxies. The ex- 
tended formula advocated here links the relative position of 
the bars to the sign of rotation of the inner bar. This exten- 
sion cannot be as rigorous as the original method, because 
it assumes that the patterns do not change as they rotate 
one through another, which is not true in general. How- 
ever, it predicts the same deviations of the integrals in the 
Tremaine- Weinberg method from their values for the single 
rotating pattern, as in the orbital model of a realistic doubly 
barred galaxy, which does not involve this assumption. This 
indicates that the degree of change that the patterns rotat- 
ing one through another undergo is small, and that it does 
not affect significantly the results of the extended method 
proposed here. 

Application of the extended method to NGC 2950 im- 
plies that the inner bar there counter-rotates with respect to 
the outer bar and to the large-scale disc. Since a retrograde 
bar is unlikely to be supported in a prograde disc, a retro- 
grade inner disc may be hiding in the central kiloparsecs of 
NGC 2950. 
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APPENDIX A: CORRECTION TERMS IN 
TREMAINE- WEINBERG INTEGRALS FOR A 
PULSATING FERRERS' BAR 

Consider a Ferrers' bar with major and minor axes a and b. 
Its surface density can be written in the Cartesian coordi- 
nates (x,y) in the plane of the galaxy as 



Y,{x,y) 







if % + 

if 4 + 



< 1 
> 1 



(Al) 



where Eos=const is the central density of the bar. Pulsa- 
tion of the bar can be described by periodic variation of the 
length of its axes 

a 2 (t) = a 2 / gi (t\P), b 2 (t) = b 2 /g 2 (t\P), 

where gi and g 2 are periodic functions of time with period 
P. Thus non-zero density, corresponding to the top line in 
(Al), of a pulsating bar that rotates with pattern speed Q,b, 
can be written in polar coordinates (R, ip) as 

E(x, y, t) = t{R, ip - Q B t, t\P) = Eos x (A2) 
[1 - MR, <p - fl B t)gi(t\P) - f 2 (R, <p - n B t)g 2 (t\P)] , 



where 

MR, 4 



MR, 



b 2 o 



Note that (A2) is one of the components of (2) for the par- 
ticular case of a Ferrers' bar. 

Let (A2), wherever larger than zero, describe non-zero 
surface density of a tracer in a rotating and pulsating bar. 
In the Tremaine- Weinberg method, the time derivative of 
this surface density enters the continuity equation. Simple 
partial derivation of (A2) gives 



ol, agi ag 2 

-dt = ~ QB d^~ EoS ^~dt~ E ° B h ~dt' 



(A3) 



Further on in the Tremaine- Weinberg method, the continu- 
ity equation is integrated over x and y, which in the case 
considered here leads to 



T,(x, y) x dx — Eos 



/i dx dy- (A4) 



Eos ^ J J f 2 dxdy = J E(x, y) v y (x, y) dx, 

where v y is the velocity of the tracer along the y axis. (A4) 
is the counterpart of (4) and (5), and since bar's pulsation 
is explicitly included here, it contains two additional cor- 
rection terms (second and third term), involving integration 
over /i and f 2 . This integration, although formally extend- 
ing to infinity, in this case is limited to the regions where 
bar's surface density is non-zero. Thus since /i and f 2 con- 
tain only rotated Cartesian coordinates, the integrals that 
involve them are finite, and will be denoted as Iib and I 2 b ■ 
Moreover, let us represent the pulsation of the bar by a sim- 
ple oscillatory form of g\ and g 2 : 



gi, 2 = 1 ± es sin tot, 



(A5) 



where lu — 2n/P, and e B controls the amplitude of the oscil- 
lation. Then, with the use of notation from Section 2, (A4) 
takes the form 

Q B L B X B - £bu cosijjt Eq B (Iib - I 2 b) = L b Vb/ sini (A6) 

Similar exercise can be done for the second bar, indexed by 
S, leading to an equation being a counterpart to (A6) 

SlsLsXs - esoj cosu)t T, s(hs ~ hs) = LsVs/ s'mi (A7) 

The sum of these two equations gives (6), the equation of the 
extended Tremaine- Weinberg method, but here with correc- 
tion terms, one for each bar. In both (A6) and (A7), the 
correction term is the second term. 
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For each bar, the integrals I\ and I2 can be evaluated 
explicitly, indicating that T,ob,s(Iib,s — l2B,s) is about twice 
smaller than Xb,sLb,s throughout each bar. Numerical sim- 
ulations indicate that the overall pulsation of the inner bar 
is more noticeable than that of the outer bar, hence the 
magnitude of the coefficient e should be larger for the inner 
bar. In the numerical N-body simulations by Rautiainen et 
al.(2002), the axial ratio of the inner bar, b/a, varies roughly 
between 0.52 and 0.72. In the notation adopted here, this 
corresponds to es = 0.16. Finally, ui = 2(Qs — Qb) is al- 
ways smaller than 2£ls- Thus the magnitude of the second, 
correction term in (A7) should not exceed about 15% of the 
magnitude of the first, leading term in the case of the in- 
ner bar. For the outer bar, the ratio of u>/Q,b is larger, but 
since the pulsation of that bar has smaller amplitude, es is 
significantly smaller. 

This example shows that periodic changes in time of 
surface density of the tracer of each bar in this bar's refer- 
ence frame can be accommodated as correction terms in the 
extended Tremaine- Weinberg formalism proposed in this pa- 
per, given that the amplitude of these changes is sufficiently 
small. This amplitude, as observed in numerical simulations, 
is indeed small enough, and it leads to correction terms that 
are one order of magnitude smaller than the leading terms 
in the extended Tremaine- Weinberg equation (6). 
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